Single-photon Optomechanics 
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Optomechanics experiments are rapidly approaching the regime where the radiation pressure of a single pho- 
ton displaces the mechanical oscillator by more than its zero-point uncertainty. We show that in this limit the 
power spectrum has multiple sidebands and that the cavity response has several resonances in the resolved- 
sideband limit. Using master-equation simulations, we also study the crossover from the weak-coupling many- 
photon to the single-photon strong-coupling regime. Finally, we find non-Gaussian steady-states of the me- 
chanical oscillator when multi-photon transitions are resonant. Our study provides the tools to detect and take 
advantage of this novel regime of optomechanics. 

PACS numbers: 42.50.Wk, 42.65.-k, 07.10.Cm, 37.30.+i 



Introduction. Optomechanics is a rapidly growing field of 
research studying mechanical degrees of freedom coupled to 
modes of optical cavities via radiation pressure, optical gradi- 
ent, or photothermal forces iH |2|. Work in this area is largely 
motivated by building more sensitive mass and force sensors 
lUt], providing long-range interaction between qubits in future 
quantum information hardware and probing quantum me- 
chanics at increasingly large mass and length scales [5]. 

In the standard optomechanics setup, the position of a me- 
chanical oscillator parametrically modulates the frequency of 
an optical cavity mode. In most experiments to date this op- 
tomechanical coupling is small compared to the mechanical 
frequency and the cavity linewidth. However, if the cavity is 
strongly driven and thus contains a large number of photons, 
the coupling between the mechanical oscillator and the fluc- 
tuations of the cavity field is enhanced by a factor y/n, where 
n is the mean photon number in the cavity. This has recently 
led to the observation of radiation-pressure effects, e^. red- 
sideband cooling di-lil, normal-mode s plit t ing 1 12, 13|, and 
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optomechanically-induced transparency ( 

In this weak coupling regime the Hamiltonian is quadratic 
so that ordinary thermal and vacuum noise lead to Gaussian 
steady-states. To create more general and possibly more inter- 
esting and useful states one either needs non-Gaussian input 
noise, e.g. driving the system with single-photon sources [|I6tl . 
or one has to make the system nonlinear. The latter can be 
achieved either via measurement backaction [17] or intrinsic 
non-linearities, e.g. coupling the resonator via a qubit to the 
mechanical oscillator [18,1 or engineeiing an optomechanical 
interaction which couples the position squared of the oscilla- 
tor to the cavity mode lloi [l9l - l2 111 . 

Several optomechanics setups, using either ultracold atoms 
in optical resonators |22], optomechanical crystals ||23| or su- 
perconducting circuits |13], are approaching the limit where 
the radiation pressure of a single photon displaces the me- 
chanical oscillator by more than its zero-point unceitainty. In 
this single-photon strong-coupling regime the full paramet- 
ric coupling, i.e. three-wave mixing in the language of quan- 
tum optics, has to be taken into account. To date there exists 
little literature on this subject with the notable exception of 
Refs. iiiH. 
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FIG. 1 : (Color online) (a) Standard optomechanics setup: the posi- 
tion a; of a mechanical oscillator is parametrically coupled to a driven 
cavity mode a. (b) Spectrum and eigenfunctions of Hamiltonian l[T). 
The energy axis is not to scale. Parabolas indicate the displaced har- 
monic oscillator potentials for n — 0,1,2 photons and g < 0. 



In this paper we show how to detect this novel regime of 
optomechanics and exploit the nonlinear spectrum to create 
non-Gaussian steady-states of the mechanical oscillator. For 
weak coherent optical drive we use the polaron transformation 
to calculate properties of the output light to all orders in the 
optomechanical coupling. We find that the power spectrum 
has multiple mechanical sidebands and the cavity response 
has several resonances in the resolved-sideband limit. Using 
master-equation simulations, we calculate these observables 
throughout the crossover from the many-photon to the single- 
photon limit. Finally, we show that multi-photon transitions 
can lead to non-Gaussian steady-states which might enable the 
observation of quantum tunneling and noise-induced switch- 
ing in optomechanical systems. 

Model. We consider the standard model of optomechanical 
systems where the position of a mechanical oscillator, x = 
xzPF{b + &^), is parametrically coupled to an optical cavity 
mode a, see Fig.[T](a). Setting h = 1 the Hamiltonian reads 



Ho = loro) a + ujMb^b + ga'a{b + b^) 



(1) 



where lor is the resonator frequency, lom the mechanical fre- 
quency, and g = oj'j^xzpf is the optomechanical coupling. 
xzPF = {2AIujm)^^^^ is the zero-point uncertainty, AI the 

mass of the mechanical oscillator, and w'„ — %^ the deiiva- 

' H ax 

tive of the resonator energy with respect to the oscillator po- 
sition X. a and b are bosonic annihilation operators for the 
cavity mode and the mechanical oscillator, respectively. 
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FIG. 2: (Color online) Detecting the single-photon strong-coupling 
regime: (a) Steady-state mean photon number (a^a) as a function of 
detuning A and (b) power spectrum S{iu) at A = 0. ljm ~ g and 
ll ~ 20 for all curves, ljm /k = 2 and nth = (blue solid), 
(jJm/i^ ~ 2 and nth ~ 1 (red dashed) as well as cjm/«: = 0.5 and 
nth = (black dash-dotted). The thin black solid line in (a) shows 
the empty cavity response g — for comparison, no is the mean 
photon number on resonance, i.e. no — AQ^ /k^ . 



Note first that the Hamiltonian ([T]i conserves photon num- 
ber, i.e. [a''a,Ho] — 0. The Hamihonian in the subspace of 
n photons is a harmonic oscillator with frequency ujm which 
is displaced by —nxo, where xq = —2xzpf9/^m is the dis- 
placement caused by one photon. Thus, the eigenvalues of ^ 
are i?„„i ~ — g^n^ I^m + i^Aim with non-negative inte- 
gers n and m. The anharmonicity is proportional to the prod- 
uct of photon number n and oscillator displacement which is 
linear in the photon number n. We show the spectrum and 
eigenfunctions of Hamiltonian ([T]i in Fig.[T](b). 

In order to include drive and decay we use standard input- 
output theory In a frame rotating at the frequency of the 
optical drive, the non-linear quantum Langevin equations read 



a = +iAa — —a — ig (^b^ -+ 
b = -iiOMb ^2^^ iga^a + V7^ii 



/«■ Sin (2) 

(3) 



where A = — lor is the detuning between laser and 
resonator frequency lur, and 7 and k are the mechanical and 
cavity damping rates. The cavity input Sin is a sum of a co- 
herent amplitude Sin and a vacuum noise operator ^ satisfy- 
ing {£.mKt')) = Sit - t') and {i^{mt')) = 0. Finally, 
we assume that the mechanical bath is Markovian and has a 
temperature T, i.e. {bin{t)bl^{t')) = (nth + l)5{t ~ t') and 
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(feL(^)&in(i')> = ntY,5{t - t') with n 

The model is characterized by three dimensionless param- 
eters: the mechanical quality factor ujmIi, the resolved- 
sideband parameter lom / k, and the granularity parameter gj n 
ll26ll28ll . The latter is the cavity frequency shift in units of its 
linewidth when the oscillator is displaced by one zero-point 
uncertainty X2P¥- Finally, 2g/ujM = '^{g/i^) x (waj/k)^^ 
is the oscillator displacement in units of xzpf caused by the 
radiation pressure of a single photon, \f\g\ > ijOm we will say 
the system is in the single-photon strong-coupling regime. 

Approximate solution for weak drive. It is well known that 
the Hamiltonian Hq can be diagonalized by the polaron trans- 
form given hyU = e'^ with S = ^StS(&t - 6) lUi]. Here 
we use it to find an approximate solution to Eqs. (|2]l and (|3]l. 



In steady-state and for weak optical drive we obtain 



a{t)^^ Sre-('^/2-^^)(*"^)e^We-*(^)Si„(T) (4) 



where we defined A = A + g"^ /lom and X{t) is given by 



X{i) 



^/lg 



i^M J- 



-iUM{t-T) 



6in(T) - H.C. 



(5) 

Using this analytic approach, we calculate properties of the 
optical field. We get for the steady-state mean photon number 
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(6) 



and the cavity amplitude relevant for homodyne experiments 



(S) 



00 

n=0 



jg/^M 

2n! 



i2n " 



E 



("th 



■th 



k=0 

(ii^) _ ,(A - (n - 2k)ujM) 



(7) 



where no = Afl^/K^ with fl = v^|Sin| is the mean pho- 
ton number for g = on resonance A = 0. The quanti- 
ties (|6]l and (|7]i are sums of resonances which are spaced by 
the mechanical frequency lum- Let us discuss first the case 
of zero temperature, i.e. nth — 0, when only terms with 
k = contribute in Eqs. ^ and (|7]i. In this case, the reso- 
nances are weighted by a Poisson distribution with variance 
{g/ujMY and the widths are n + n^. The resonances can be 
understood as transitions between the vacuum state |0, 0) and 
the manifold of one-photon eigenstates |1, m) of the Hamilto- 
nian ([U. They are resonant if the laser frequency cjl matches 
bJL = Eim^Eoo = ojR — g'^ /ujm + 'mujM . The Poission dis- 
tribution is due to the Franck-Condon factors Kmje'^lO)!^ = 
1/ dxip;„{x~XQ)ipoix)\^ = (.g/wA/)2"e-(9/"'>^)V"^! where 
|m) is the state with m phonons, (pm (x) is its real-space wave- 
function, and xo/xzPF = 2g/uJM- At finite temperature the 
states |0, m) with m > are thermally occupied leading to a 
redistribution of weight among the peaks and additional reso- 
nances at ujl — ujji — g^ juu — mujM- 

In the limit k ^ 7 we obtain the cavity spectrum S{uj) — 
J^^dte^-^[{aHt)aiO))-\{a{tm as 



Cmn{m + n)7 



771,71 — 



(8) 



with the n = rn = term excluded. The coefficients C,„„ are 
independent of cj but rather involved and will not be shown. 

The optical output spectrum has sideband peaks at integer 
multiples of the mechanical frequency lum whose widths are 
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FIG. 3: (Color online) Crossover from many- to single photon limit: 
Steady-state mean photon number (a^a) as a function of detuning A. 
(a) Q./k = 0.01 and q/k = 2, (b) Q/k ^ 0.5 and g/K. = 0.5, and 
(c) Q/k — 20 and g/n = 0.01. Parameters are ujm/h — 2, A/k = 
—2, ljai/j = 100 and nth = 0. We show Eq. JSJ (red solid), jap 
(blue dash-dotted), linear theory |ap + (d^d) (green dashed) and 
simulations of Eq. ^ (black dots), (d-f) Output spectrum S'(cj) at 
A = for the same parameters from simulations of Eq. (|9]. 



multiples of the mechanical linewidth 7. At zero temperature 
there are peaks only at negative frequencies because photons 
can only create phonons and leave the cavity with frequencies 
smaller than the laser frequency lul. At finite temperature 
(or stronger optical drive) additional peaks appear at positive 
frequencies since there is a finite probability that a photon ab- 
sorbs the energy of one or more phonons and leaves the cavity 
with a frequency larger than the laser frequency w^. In pass- 
ing we note that driving on these additional sidebands in the 
resolved-sideband limit leads to multi-phonon cooling which 
will be discussed in a future publication. 

In Fig. |2] (a) we plot the steady-state mean photon num- 
ber (a^a) as a function of detuning A for a system entering 
the single-photon strong-coupling regime, g — lum. In the 
good-cavity limit k < ljm the cavity response shows several 
resolved resonances. At finite thermal phonon number addi- 
tional peaks appear and their weights are redistributed until 
eventually they blur into a broad thermal background. In the 
bad-cavity limit k > ljm the resonances overlap and broaden 
the empty cavity resonance. In Fig.|2](b) we present the out- 
put spectrum S{u!) at zero detuning A = 0. It shows a series 
of peaks at multiples of the mechanical frequency ojm for all 
sets of parameters: within and outside the resolved-sideband 
limit as well as for a zero and finite thermal phonon number 

Crossover between the many- and the single-photon limit. 
Let us now compare the single-photon strong-coupling regime 
to the more familiar case of weak optomechanical coupling 
and study the crossover between these two extreme limits. 

For numerical simulations it is advantageous to use the dis- 
placement transformation by writing a = a + d and b = 
6 + c. We obtain a set of coupled equations for the mean 



values a and h: 



zAa 



i^l — ig{h + h*)a and 



= —iijjMh—^h — ig\a\ . It is well-known that these nonlin- 
ear equations have either one or three solutions. In the latter 
case the system is said to be (classically) bistable. The opera- 



tors d and c describing the fluctuations around the mean values 
a and b, respectively, satisfy equations of motion equivalent to 
the quantum master equation 

H'^,Q + KV[d]g + "f[nt\-, + l)V[c]Q + -^nt\J)[c^]Q 

(9) 

with the Hamiltonian 

//q = - /\' d+LOMc) c+g{d* d+dd)) (c + c^)+gd'^d (c + c^) 

(10) 

where the detuning is renormalized A' — /S. — g{h + h*). 
T)[o\q = OQO^ — {o^dg + go^d)/2 is the standard dissipator in 
Lindblad form. This is an exact description of the system in a 
frame where the mean of both harmonic oscillators has been 
displaced to the vacuum. 

Outside the bistable region, for large mean cavity amplitude 
d and small optomechanical coupling g, the last term in the 
Hamiltonian ( fTOb can be neglected. We then obtain a quadratic 
Hamiltonian or equivalently a set of linear quantum Langevin 
equations which can be solved exactly. In this linear theory we 
have {d) —{), and the photon number is given by | a p + (d^ d) . 

We now compare the predictions of the numerical solu- 
tion of the quantum master equation (|9]l to the linear the- 
ory and the analytic expressions (|6) and (|8) derived above. 
In Fig. [3] (a-c) we plot the steady-state mean photon num- 
ber (a^d) as a function of detuning A for three different sets 
of parameters. For Vl/n = 0.01 and g/n = 2 we are in 
the single-photon strong-coupling limit. The numerical so- 
lution of Eq. (|9]l shows several resonances and agrees very 
well with the analytical expression (|6]l. The linear theory is 
not appropriate in this regime. This is signaled by the fact that 
the size of the fluctuations by far exceeds the mean photon 
number: |ap ^ {d^d). At intermediate coupling and drive, 
Q,/ K = 0.5 and g/ k = 0.5, the numerical simulation of ^ 
predicts one large peak slightly below A = and a small 
resonance close to the blue sideband A = lom- Eq. ^ and 
the linear theory qualitatively describe this feature but fail to 
agree with the numerics quantitatively. Finally, for VL/ k, — 2Q 
and g/ K — 0.01 we are well inside the regime where the lin- 
ear theory is vaUd. It correctly predicts a slightly asymmetric 
peak close to A = 0. 

In Fig.[3](d-f) we show the optical output spectra 5(a;) for 
the same parameters obtained from simulations of Eq. (|9}. In 
the single-photon strong-coupling limit it has multiple side- 
bands and agrees quantitatively with Eq. (|8]l. As the drive 
strength increases, additional sidebands at positive frequen- 
cies appear. With decreasing optomechanical coupling g the 
weight gradually concentrates in the two sidebands at = 
±ijOm as predicted by the linear theory. 

Non-Gaussian steady-states via multi-photon transitions. 
In the final part of this paper we address the question as to 
how the nonlinear single-photon strong-coupling limit leads 
to non-Gaussian steady-states for the mechanical oscillator. 

Recalling the spectrum of the Hamiltonian ([T]) we notice 
that for A = —ng'^/uJM multi -photon transitions between the 
vacuum state |0, 0) and the lowest-energy state with n pho- 
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FIG. 4: (Color online) Non-Gaussian steady-states via multi-photon 
transitions, (a) Steady-state mean phonon number (b^b) (blue solid) 
and the second-order coherence of the mechanical oscillator F — 
{b''b''bb) /{{b'b))^ (green dashed) as a function of drive strength Q,. 
(b) Phonon number distribution Pn at SI/k — 0.6. Parameters are 



A : 



-3g Ilom,^m 1 ~ 2,ujM/'y ~ 1000, andg/K = 1. 



tons |ri, 0) are resonant, i.e. Eqq + ulol — Eno- As all in- 
termediate transitions are off-resonant, the system is in the 
Franck-Condon blockade regime fs^, and we expect that at 
weak drive the system will stay close to the vacuum state and 
a strong drive induces multi-photon transitions. 

In Fig.|4](a) we show the steady-state mean phonon num- 
ber n = (6^6) and the second-order coherence of the mechan- 
ical oscillator F = {b^ w bb) / {{b'b))'^ as a function of drive 
strength fl. For small drive rt/n < 0.3 the mean phonon 
number n remains small as expected. After a region of strong 
number fluctuations F ^ 1, the phonon number increases 
rapidly with drive strength il. We note that the phonon num- 
ber fluctuations also become large in the limit of small drive 

— > which has been discussed in Ref. Isoll . 

Since | (6) | is small for all drive strengths ft considered here 
(not shown), the Wigner function of the mechanical oscillator 
is rotationally invariant and the phonon number distribution 
Pn contains the complete information of the reduced density 
matrix. In Fig. |4](b) we plot the phonon number distribution 
Pn for n/n = 0.6. We can clearly distinguish two peaks: one 
at zero and one at n = 14 phonons. With increasing drive 
strength il weight is gradually transferred from the former to 
the latter We interpret this as a statistical mixture of two dif- 
ferent oscillation amplitudes and expect the system to exhibit 
quantum tunneling and noise-induced switching [31]. We note 
that similar states have been reported in Ref. [26]. 

Conclusion. Motivated by recent experiments, we explored 
optomechanics in the regime where the radiation pressure of 
a single photon displaces the mechanical oscillator by more 
than its zero-point uncertainty. We demonstrated that the out- 
put spectrum and cavity response are qualitatively modified 
and showed how to create non-Gaussian steady states of the 
mechanical oscillator. Our study opens many further ques- 
tions about the single-photon strong-coupling regime whose 
physics is far from well-understood, including e.g. the fate of 
red-sideband cooling and ponderomotive squeezing. 
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Note added. During the final stages of this project, a related 
paper by Rabl appeared ||32 ] . 
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